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Abstract. We consider the real part Re C(s) of the Riemann zeta- 
function C(s) m the half-plane Re (s) > 1. We show how to com- 
pute accurately the constant ^ 1.19 which is defined to be the 
supremum of tr such that Re C,{u + it) can be negative (or zero) for 
some real t. We also consider intervals where ReC(l + it) <Q and 
show that they are rare. The first occurs for t k, 682112.9, and has 
length « 0.05. We list the first 50 such intervals. 



1. Introduction 

In this note we consider the real part of the Riemann zeta-function 
in the half-plane if = {s G C | Re (s) > 1}. As usual, we write 
s = a + it, so Re (s) = cr > 1. We are mainly interested in the 
regions where ReC(s) < 0. Since limo-t-oo C(^ + ^^) = 1 (uniformly in t), 
ReC(o" + it) cannot be zero for arbitrarily large a > 1. We define 

(To := sup{(T e M I (3t G M) Re C((T + it) = 0}. 

Thus, ReC(s) > if cr > ctq. In van de Lune [9] it was shown that 
(Jo is the (unique) positive real root of the equation 

Earcsinf — ^ = — , 

p 

where p runs through the primes (we adopt this convention through- 
out). In [9] it was also shown that ctq > 1.192 and that Re(^(o"o + it) 
never vanishes. 

The main aim of this note is to show how ctq can be computed to 
arbitrarily high precision by an efficient algorithm. We also mention 
some results on the behaviour of ReC(cr -|- it) for 1 < a < aQ, and in 
particular on the line a = 1. 

2. Accurate computation of the constant ctq 

In this section we assume that a > ai > 1, where ai is a suitable 
constant (e.g. 1.1). We show how the constant (Tq can be computed 
within a given error bound. There are three main steps. 

1 
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(1) Give an algorithm to evaluate the prime zeta-function [5] 

p 

for real a > 1. 

(2) Using step 1, give an algorithm to evaluate the function /(a) 
defined by 

/(a) = ^arcsin(^) 

p 

(3) Use a suitable zero-finding algorithm to locate a zero of f{a) 
in a (sufficiently small) interval where f{a) changes sign, for 
example [1.1, 1.2]. 

Step 1 is easy. From the Euler product for ({a) and Mobius inversion, 
we have a formula essentially known to Euler [4, 1748]: 

(1) P(a) = f;^logCM, 

r=l 

which is valid for cr > 1 (see Titchmarsh [13, eqn. (1.6.1)]). The series 
converges rapidly in view of the following Lemma. 

Lemma 2.1. Fora>2,0 < logC(a) < 3/2'^ and < P{a) < 3/2'^. 

Proof. For cr > 2, we have 

poo Ol — (T 

< C{a) - 1 < 2"'^ + 3-^" + / x-'^dx = 2"^^ + 3"^" + < 3/2'", 

J 3 cr - 1 

SO 

0<logC(o^)<C(^)-l<3/2-. 
The upper bound on P{cr) follows similarly, using P(cr) < ([a) — 1. □ 

Using (1) and Lemma 2.1, we have 

/i(r) 



oo 



P{a) = logC(a) + 2^ ^ logC(ra) 

r=2 ^ 



where the r-th term in the sum is bounded in absolute value by 3/2'"'"+^. 
Thus, we can evaluate P{(j) accurately, for given cr > 1, using any 
good algorithm for the evaluation of C('^)) ioi example Euler-Maclaurin 
summation. If (1) is used to compute -P(cr), P(3cr), P(5cr), . . . , then 
we should take care to compute the relevant terms log ({ra) only once. 

For step 2, we observe that the arcsin series defining f{a) converges 
slowly and irregularly, since it is a sum over primes which to first or- 
der behaves like YlpP"'^- The well-known "trick" is to express /(cr) 
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as a double series and reverse the order of summation, obtaining an 
expression which is mathematically equivalent but computationally far 
superior. For some similar examples, see Wrench [15, 1961]. 
For |x| < 1 we have 

oo 

arcsin(x) = Cfca:^^"'"\ 

fc=0 

where 

l-3-5---(2A;-l) 1 (2A;)! , , 

2-4-6---(2A;) 2A; + 1 (2'=A;!)2(2fc + 1) «^ - ^• 

Note that all Cfc are positive so that /(cr) is strictly convex. It is also 
clear that /(cr) is strictly decreasing for c > 1. From the expression 
for Cfc, we see that, for > 1, 

For cr > 1 it is easy to justify interchanging the order of summation 

in 

oo ^ 
p fc=0 

obtaining 



1 \2fc+l TT 

.p^/ ~ 2 



(3) / w = E E - i = E ^^^((2^ + 1)-) - \ ■ 

k=0 p ^ k=0 

From Lemma 2.1 and the inequality (2), we see that 

oo 

0< E CfcP((2A; + l)a) <2-(2^+3)<x^ 

k=K+l 

SO it is easy to determine K such that we can truncate the series in (3) 
to a finite sum over k < K with a rigorous error bound. 

If desired, we can substitute (1) into (3) and interchange the order 
of summation, obtaining 

(4) /(a) = E^.logC((2j + lV)-^, 

i>o 

where 

Cfc/u(r) 



r 

fc>0,r>l,(2/c+l)r=2j+l 
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From the inequality Ck < l/{2k + 1) (valid for k > 0), it follows that 
\dj\ < 1- Using Lemma 2.1, we can determine where to safely truncate 
the series (4). 

For step 3, we can use a zero-finding algorithm which needs only func- 
tion (not derivative) evaluations, and gives a guaranteed bound on the 
final result. For example, the method of bisection could be used, but 
would be slow, taking about log2(l/£:) function evaluations to obtain 
a solution with error bounded by e. In the secant method, a sequence 
{xn), converging to a zero of / under suitable conditions, is obtained 
by computing the approximation x„ by linear interpolation using the 
two points (x„_i, /(x„„i)) and (x„_2, f{xn-2))- It converges with order 
(1 + -\/5)/2 ~ 1.618, but does not always give a guaranteed bound on 
the error. A combination of bisection and linear interpolation, as in 
the algorithms of Dekker [3] or Brent [2], can give convergence about 
as fast as the secant method, but with the final result bracketed in a 
short interval where the function / changes sign. 

3. Computational results 

The second and third authors independently wrote programs imple- 
menting the ideas of §2, using Magma in one case and Mathematica 4 
and 8 in the other case. The programs used different strategies to ob- 
tain a final interval where / changes sign (in one case taking advantage 
of the strict convexity of /). The output of the programs agrees to at 
least 500D. We give here the correctly rounded result to lOOD: 



ao ^ 1.19234 7337186193 20289 7504427425 59788 34011 19230 83799 
94301 37194 92990 52458 64848 30139 24084 99863 83788 36244 . 

Programs and higher precision values are available from the authors. 

4. The distribution of ReC(cr + it) for a > 1 
Assuming that the limit exists, we define 

d{a)= lim m{t G [0, T] | ReC(cT + it) < 0}, 

where m denotes Lebesgue measure. Informally, d{a) is the probability 
that C(s) has negative real part on a given vertical line Re (s) = a. 

The results of Section 2 show that d{a) = for cr > (Tq ~ 1.19. Here 
we briefly discuss the region 1 < a < (Jq. 

At least for those values of t that are accessible to computation, 
ReC((j -|- it) is "usually" positive for a > 1. The function d{a) is 
conjectured to be continuous and monotonic decreasing from a positive 
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value at cr = 1 to zero at a = ctq- Even on the line a = 1, Re C (a + it) is 
usually positive [11]. We can prove that d{l) < 1/4, but a Monte Carlo 
computation suggests that the true value is much smaller. Based on 
5 X 10^^ pseudo-random trials, we estimate d{l) = (3.80 ±0.01) x 10"''. 
Similarly, we estimate rf(l.Ol) = (1.10 ± 0.01) x 10"^ and c/(1.02) ^ 
(2.66 ± 0.04) X 10~^, so it can be seen that d{a) decreases rapidly as 
we move to the right of a = 1. 

Although ({s) has a simple pole at s = 1, the Laurent series 

C(5) = ^ + 7 + 0(|s-l|) 
s — 1 

shows that Re(^(l + it) has a positive limit 7 = 0.577- ■■ (Euler's 
constant) as t — )■ 0. 

On any fixed vertical line a > 1, both ({a + it) and l/C(o" + "it) 
are bounded, in fact ({2a)/({(r) < \({(t + it)\ < ({a). However, the 
situation is different on the line a = 1, as both ({1 + it) and + it) 
are unbounded. Their true order of growth is unknown. It follows from 
Titchmarsh [13, Theorem 11.9] and the continuity of ReC(l + it) that 
lie ({1 + it) attains all real values. Nevertheless, the "usual" values are 
quite small. As a special case of [13, Theorem 7.2] we have the mean 
value theorem 

hm i ^\al + ^t)\'dt = a2) = ^. 

T-s-oo I Ji 6 

Using ideas as in the proof of [13, Theorem 7.2], we can prove that 

1 /"^ 

lim - / Re((l + it) dt = 1. 

Thus, informally, we can say that the typical value of Re(^(l + it) is 
close to 1. The values have a distribution with mean 1 and variance 
7rV6 - 1 ^ 0.645. 

In [9, Table 1], van de Lune gave a list of values of t > such 
that ReC(l + it) < and is (approximately) a local minimum. The 
list was not claimed to be exhaustive. The smallest t listed was t = 
682112.92 with ReC(l + it) ^ -0.003. We have shown, using the 
"maximum slope principle" [10], that this is very close to the smallest 
t for which ReC(l + it) < 0. More precisely, ReC(l + it) > for 
< t < 682112.8913, and there is a local minimum of -.0027652 at 
t ^ 682112.9169. In applying the maximum slope principle we used 
the bound 

< ^ log{t^ + 4) + 7 for t> 10. 



dt 



argC(l + it) 



Re 



C(l + ^t) 

C(l + It) 
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Table 1. First 50 negative local minima of ReC(l + it) 



t 


ReC 


length 


t 


ReC 


length 


682112.9169 


-0 


0028 


0.0529 


8350473.4853 


-0 


0019 





0451 


1267065.1710 


-0 


0040 


0.0655 


8366684.0439 


-0 


0197 





1322 


1466782.0667 


-0 


0013 


0.0391 


8452317.9526 


-0 


0090 





0900 


1858650.0915 


-0 


0282 


0.1686 


8967566.5926 


-0 


0148 





1336 


2023654.7671 


-0 


0221 


0.1389 


9960968.8748 


-0 


0184 





1373 


2064996.2141 


-0 


0117 


0.1076 


11231380.7309 


-0 


0099 





1042 


2195056.7909 


-0 


0755 


0.2718 


11236680.3350 


-0 


0262 





1595 


2202620.3296 


-0 


0111 


0.1159 


11781932.0257 


-0 


0170 





1288 


2530662.6360 


-0 


0072 


0.0865 


11884021.9776 


-0 


0035 





0564 


3259774.5293 


-0 


0471 


0.2098 


12045289.3337 


-0 


0644 





2498 


3548283.4160 


-0 


0189 


0.1459 


12276788.1573 


-0 


0182 





1476 


4052438.9330 


-0 


0023 


0.0474 


12546625.7916 


-0 


0455 





2031 


4197235.0783 


-0 


0331 


0.1977 


12781127.5748 


-0 


0102 





0964 


5410820.7150 


-0 


0008 


0.0307 


13598773.5889 


-0 


0543 





2317 


6027913.8513 


-0 


0181 


0.1325 


13786262.5457 


-0 


0826 





2635 


6164063.0008 


-0 


0263 


0.1603 


13922411.7750 


-0 


0222 





1418 


6238849.4877 


-0 


0071 


0.0827 


14190358.4974 


-0 


0632 





2214 


6265907.4688 


-0 


0030 


0.0522 


14391623.0217 


-0 


0016 





0437 


6421627.2235 


-0 


0241 


0.1651 


14788310.5330 


-0 


0149 





1132 


7338152.4379 


-0 


0043 


0.0656 


14856540.3430 


-0 


0220 





1442 


7469838.9709 


-0 


0009 


0.0305 


15173904.7533 


-0 


0041 





0800 


7766995.0303 


-0 


0742 


0.2840 


15321273.7219 


-0 


0131 





1181 


7774558.3985 


-0 


0672 


0.2705 


16083163.0244 


-0 


0098 





1038 


7985493.9836 


-0 


0324 


0.1728 


16503899.3235 


-0 


0060 





0680 


8299958.2327 


-0 


0022 


0.0432 


16656258.8346 


-0 


0155 





1329 



Table 1 lists the first 50 local minima of ReC(l + it) for which 
t > and ReC(l + it) < (no minima are exactly zero). The val- 
ues in the table are rounded to 4 decimal places. The columns headed 
"length" give the lengths of the intervals containing t in which Re C is 
negative. To 8 decimal places, the first interval, of length 0.05291225, 
is (682112.89133824,682112.94425049). The sum of the lengths of the 
first 50 intervals is 6.48390168, giving an estimate d{l) ~ 3.85 x 10"'^. 
This is close to our Monte Carlo estimate d{l) ~ 3.80 x 10~^. 

In this brief note we refrain from commenting on the region 
a G [1/2,1), but refer the interested reader to the literature, such 
as Bohr and Jessen [1], Titchmarsh [13, §11.13], Tsang [14], Joyner [6], 
Laurincikas [8], Steuding [12] and Kiihn [7]. 
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